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Using the thermal quasiparticle random phase approximation, we study the process of neutrino 
and antineutrino capture on hot nuclei in supernova environments. For the sample nuclei 06 Fe 
and 82 Ge we perform a detailed analysis of thermal effects on the strength distribution of allowed 
Gamow-Teller transitions which dominate low-energy charged-current neutrino reactions. The finite 
temperature cross sections are calculated taking into account the contributions of both allowed and 
forbidden transitions. The enhancement of the low-energy cross sections is explained by considering 
thermal effects on the GT± strength. For 56 Fe we compare the calculated finite-temperature cross 
sections with those obtained from large-scale shell-model calculations. 
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I. INTRODUCTION 

It is well known that core-collapse supernova simula¬ 
tions require a detailed description of neutrino transport 
including all potentially important neutrino reactions [1], 
Neutrinos are the mediators of the energy transfer from 
the core to the outer stellar layers and their luminosities 
and spectra are a crucial ingredient for the supernova 
explosion mechanism. Despite significant progress in our 
understanding of the core-collapse mechanism, many su¬ 
pernova simulations fail to produce explosion (see e.g. 
[2, 3] and references therein). One of the possible reasons 
for this failure could be the incomplete and inaccurate 
treatment of the neutrino-nucleus processes occurring in 
supernova environments. 

It was pointed by Haxton [4] that neutral- and charged- 
current neutrino reactions on nuclei involving the excita¬ 
tion of the giant resonances might help reenergize the 
explosion and they should, therefore, be taken into ac¬ 
count in supernova simulations. Such simulations were 
performed by Bruenn and Haxton [5] and it was shown 
that inelastic neutrino-nucleus scattering plays an impor¬ 
tant role in equilibrating neutrinos with matter. How¬ 
ever, the simulations did not confirm Haxton’s sugges¬ 
tion that neutrino-nucleus reactions might help to revive 
the stalled shock wave. 

In [5], the nuclear composition of the core was approx¬ 
imated by a single representative nucleus, 56 Fe, and the 
relevant cross sections and rates were evaluated consider¬ 
ing only allowed (Gamow-Teller) and first-forbidden up¬ 
ward transitions from the nuclear ground state. How¬ 
ever, in the hot supernova environment with tempera¬ 
tures T > 1 MeV nuclei exist as a thermal ensemble, i.e., 
nuclear excited states are thermally populated according 
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to the Boltzmann distribution. As was first realized in 
Ref. [6], downward transitions from excited states remove 
the reaction threshold and can significantly enhance the 
reaction cross section at low neutrino energies. 

From a microscopic point of view, there are two ways 
how to treat neutrino reactions with hot nuclei. The first 
one involves a state-by-state summation over Boltzmann- 
weighted, individually determined contributions from nu¬ 
clear ground and excited states. To apply this method, 
one needs to know the strength distribution of electro- 
weak transition operators for thermally populated nu¬ 
clear states. The second method is based on a statisti¬ 
cal formulation of the nuclear many-body problem. In 
this approach rather than computing individual strength 
distributions one determines an "average" temperature- 
dependent strength function. 

For iron group nuclei (A = 45 — 65) the first approach 
was developed in [7 9] on the basis of large-scale shell- 
model (LSSM) diagonalization calculations. Modern 
high-performance compute capabilities combined with 
state-of-the-art diagonalization approaches are able to 
provide detailed strength distributions for both charge- 
neutral and charge-changing Gamow-Teller transitions 
that dominate neutrino reactions with p/-shell iron group 
nuclei at low neutrino energies {E v < 20 MeV). However, 
at temperatures T > 1 MeV an explicit state-by-state 
summation over all thermally populated states presently 
remains computationally prohibitive. To overcome this 
problem the "Brink hypothesis" is applied, i.e., it is as¬ 
sumed that GT distributions on nuclear excited states are 
the same as for the ground state. To account for ther¬ 
mal effects, the giant resonances built on the ground and 
low-lying daughter states in the inverse reaction are con¬ 
sidered. Such excited states are called back-resonances 
and their importance arises from the large nuclear ma¬ 
trix elements and increased phase space. Thus, within 
LSSM calculations the temperature dependence of the 
cross sections is comprised solely in the back-resonance 
terms. 

Within the shell model approach it is assumed that 
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the GT distributions on thermally excited states are 
the same as for the ground state. This is not likely 
the case as the vanishing of pairing correlations and 
smearing of the Fermi surface with increasing temper¬ 
ature should affect the distribution centroid and move it 
slightly down in energy. This conjecture is confirmed by 
shell model Monte Carlo studies performed at finite tem¬ 
perature [10]. In addition, with present compute capa¬ 
bilities the shell-model diagonalization cannot be applied 
for nuclei beyond the p/-shell because of the huge con¬ 
figurational spacse involved. These shortcomings can be 
avoided in a thermal quasiparticle random-phase approx¬ 
imation (TQRPA). In [11], the TQRPA was proposed as 
a method to study the response of hot nuclei to weak ex¬ 
ternal perturbations in the framework of a statistical ap¬ 
proach. Being based on the thermo-field dynamics (TFD) 
formalism [12-14], the TQRPA enables the computa¬ 
tion of temperature-dependent strength function avoid¬ 
ing the assumption of Brink’s hypothesis. In [15, 16], 
the TQRPA was applied to study thermal effects on 
neutral-current inelastic neutrino scattering on 54 ’ 56 Fe 
and neutron-rich germanium isotope 82 Ge. It was shown 
that the TQRPA reveals the same thermal effects on the 
cross sections as the LSSM approach. Namely, the re¬ 
action threshold for inelastic neutrino-nucleus scattering 
is removed at finite temperature and the cross section 
for low-energy neutrinos is significantly enhanced. It was 
found, however, that within the TQRPA the enhance¬ 
ment is caused by both de-excitation of nuclear excited 
states and thermally unblocked low-energy GT transi¬ 
tions. The latter do not appear within the LSSM due 
to application of Brink’s hypothesis. Moreover, it was 
shown in [15] that, unlike in the LSSM approach, the 
principle of detailed balance is not violated within the 
TQRPA and it results in a larger strength for downward 
GT transitions from excited states. As a consequence, 
at low neutrino energies the finite temperature cross sec¬ 
tions calculated within the TQRPA turn out to be several 
times larger than those obtained within the shell-model 
calculations. 

In the present paper, we apply the TQRPA method 
to study thermal effects on charged-current neutrino- 
nucleus reactions ( v e ,e~ ) and (V e ,e + ) occurring in the 
supernova environment. Here we would like to men¬ 
tion that the TQRPA was already applied for charge¬ 
changing reactions with hot nuclei when studying ther¬ 
mal effects on electron capture in supernovae [17, 18]. 
The paper is organized as follows: in Sec. II we review 
the basics of the TFD formalism and outline how to treat 
charge-changing transitions in a hot nucleus within the 
TQRPA. For more details of the approach, the reader 
is referred to Refs. [17, 18]. In addition, in Sec. II the 
expressions necessary to calculate cross-sections for v e - 
and 77 e -absorption on hot nuclei are given. The results of 
the numerical calculations are presented and discussed 
in Sec. Ill for the sample nuclei 56 Fe and 82 Ge. For 
56 Fe we compare the calculated ground-state and finite- 
temperature cross sections with available results from 


other approaches. Conclusions are drawn in Sec. IV. 
In an Appendix we prove in a model independent way 
that the principle of detailed balance is valid for charge¬ 
changing transitions if hot nuclei in the supernova envi¬ 
ronment are treated in the grand canonical ensemble. 

II. FORMALISM 

During the core-collapse phase of a supernovae ex¬ 
plosion the temperature in the iron core is sufficiently 
high (a few 10 9 K) to establish an equilibrium of reac¬ 
tions mediated by the strong and electro-weak interac¬ 
tion [1]. Neglecting weak-interaction reactions, we can 
consider nuclei as open quantum systems in thermal equi¬ 
librium with heat and particle reservoirs and, hence, they 
can be described as a thermal grand canonical ensem¬ 
ble with temperature T and proton and neutron chem¬ 
ical potentials X p and X n , respectively. In TFD, such 
ensemble is represented by the thermal vacuum |0(T)), 
which is a temperature-dependent state in the extended 
Hilbert space [19]. The thermal vacuum is determined as 
the zero-energy eigenstate of the thermal Hamiltonian, 
7~L = H — H , and it satisfies the thermal state condition 

A\0(T)) = a A e' H/2T Ai\0(T)). (1) 

In the above equations H is the original nuclear Hamil¬ 
tonian (proton and neutron chemical potentials are in¬ 
cluded in H) and H is its tilde counterpart acting in 
the auxiliary Hilbert space; an operator A acts in the 
physical Hilbert space, A is its tilde partner, and a a is 
a phase factor. The thermal state condition guarantees 
that the expectation value (0(T)|A|0(T)) is equal to the 
grand canonical average of A. In this sense, relation (1) 
is equivalent to the Kubo-Martin-Schwinger condition for 
an equilibrium grand canonical density matrix [20]. 

External perturbations mediated by weak interaction 
induce transitions from the thermal vacuum to nonequi¬ 
librium states. Within the TDF such nonequilibrium 
states are given by the eigenstates of the thermal Hamil¬ 
tonian "H. As follows from the definition of "H, each of its 
eigenstates with positive energy has a counterpart — the 
tilde-conjugate eigenstate — with negative but the same 
absolute value of energy. This gives the possibility to 
describe both endoergic and exoergic neutrino reactions 
with hot nuclei. In what follows we will refer to positive 
energy eigenstates as non-tilde states, and to negative 
energy eigenstates as tilde states. 

A. Thermal quasiparticle RPA 

In [17, 18] we have introduced the proton-neutron 
TQRPA method which allows for a treatment of 
charge-changing transitions in hot nuclei induced by 
(anti)neutrino absorption. For the sake of completeness, 
let us briefly recall the method. 
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Within the TQRPA, nonequilibrium states of a hot 
nucleus caused by an external perturbation are treated 
as phonon-like excitations on the thermal vacuum 

\QjMi) = 

\QjMi)=Qy.\0(T)), (2) 

where we denote = (—1) J ~ M Q\_ Mi - Phonon exci¬ 
tations are considered as the normal modes of the thermal 
Hamiltonian: 

T-L ~ ^ uji(T)(Q\ Mi Q JMi — Q^jMiQjMi)i ( 3 ) 

J Mi 

while the thermal vacuum |0(T)) itself is the vacuum for 
the annihilation operators Qjmi-, Qjmi • Thus, within 
the TQRPA the problem of finding the excitation spec¬ 
trum of a hot nucleus is reduced to the diagonalization 
of the thermal Hamiltonian in terms of phonon opera¬ 
tors such that the respective phonon vacuum obeys the 
thermal state condition (1). 

For charge-changing multipole transitions the phonon 
operators are defined as a linear superposition of cre¬ 
ation and annihilation operators of proton-neutron ther¬ 
mal quasiparticle pairs [21] 

QjMi=^2{ i’jpjn Wj p Pjn} M + i’jlip M 

jpjn 



In turn, thermal quasiparticles are normal modes of the 
pairing part of the thermal Hamiltonian 

^pair — e o(T)(PjmPjm — Pj m Pjm), (5) 

r jm 

and their vacuum is the thermal vacuum in the BCS ap¬ 
proximation. In the expression above, the notation 
implies a summation over neutron (t = n) or proton 
(r = p) single particle states only. Thermal quasipar¬ 
ticles are connected with Bogoliubov quasiparticles by 
so-called thermal transformation [22] 

ftjm. ~ x j a jm ~ 

Pjm = X 0^]m + iyjVjrn, ( x j + Vj = !)• (6) 

The (x, ^[-coefficients as well as the (it, ^[-coefficients of 
the Bogoliubov transformations are found from the finite- 
temperature BCS equations (see [18] for more details). In 
accordance with the BCS theory [23, 24], the numerical 
solution of these equations yields vanishing pairing cor¬ 
relations above a certain critical temperature. 

To clarify the physical meaning of different terms 
in (4), we note that the creation of a negative energy 


tilde thermal quasiparticle corresponds to the annihila¬ 
tion of a thermally excited Bogoliubov quasiparticle or, 
which is the same, to the creation of a quasihole state 
(see [18] for more details). Therefore, at finite temper¬ 
ature, single-particle charge-changing transitions involve 
excitations of three types: 1) two-quasiparticle excita¬ 
tions described by the operator /?t Bj n and having energy 
= Sj p + £j n , 2) one-quasiparticle one-quasihole ex¬ 
citations described by the operators B>\ , B\ B\ and 

° Jp Jn ' Jp 1 Jn 

having energies £j~j n = £j p — £j n an d — £ j~j n i respectively, 
and 3) two-quasihole excitations described by the opera¬ 
tor P] p P] n and having energy ~£j p j n - The last two types 
are possible only at T ^ 0. Therefore, due to single¬ 
particle transitions involving annihilation of thermally 
excited Bogoliubov quasiparticles, the phonon spectrum 
at finite temperature contains negative- and low-energy 
states which do not exist at zero temperature and these 
"new" phonon states can be interpreted as thermally un¬ 
blocked transitions between nuclear excited states. 

To find the energy and the wavefunctions of thermal 
phonons we apply the equation of motion method 

(\SQ,[H,Q']]\)=u(T)(\[6Q,Qi}\) (7) 

under two additional constraints: (a) the phonon vacuum 
obeys the thermal state condition (1), and (b) phonon 
operators obey Bose commutation relations. The first 
constraint yields the following relations between phonon 
amplitudes 

Jl = yjpVu - e ~ ulJil ‘ 2Tx j P x j P (<t>\ Jl 

\0/ JpJn e-“ Ji/2T yj P yi n - x j P x 3n \Wj p j n ’ 

(v\ Jl _ yj P x j n — e~ UJil2T Xj p yj n / A Jl 
\dj r j n e- UJi/ 2 T yj P x Jn-Xj p y jn \v) jpj J 

while the last assumption is equivalent to averaging with 
respect to the BCS thermal vacuum in the equations 
of motion (7) and leads to an orthonormality condition 
for the amplitudes [18]. Thus, the phonon energies oj 
in Eq.(3) as well as the amplitudes ip, ip, etc. in Eq.(4) 
are the solution of the proton-neutron TQRPA equa¬ 
tions. Positive energy solutions correspond to non-tilde 
phonons in (3), while negative energy solutions corre¬ 
spond to tilde ones. Since both the energies of ther¬ 
mal quasiparticles and the interaction strengths between 
them are temperature dependent, the spectrum of ther¬ 
mal phonons turns out to be temperature dependent. 
However, it is significant that in the zero-temperature 
limit the described method reduces into the standard 
proton-neutron QRPA. 

For a given multipole charge-changing transition op¬ 
erator T}^ the transition probabilities (strengths) from 
the thermal vacuum to thermal one-phonon states are 
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given by the following reduced matrix elements 

^ = \(Qji\\Tj T) mT))\ 2 , 

^ = \(Qji\\rj T) \m)\ 2 . 


(9) 


vacuum to a thermal one-phonon state is given by 
daji{E v ,T) ( G f cos9 c ) 2 


(El 
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xF(±Z+l,E e )(l-f(E e )) 


(13) 


Here the symbol (—) refers to neutron-to-proton transi¬ 
tions (/?“ channel), while (+) refers to proton-to-neutron 
transitions (/3 + channel). It can be shown that if the 
transition operators Tj~^ and 7j + ^ differ only by the 
isospin operator, i.e., Tj^ = 7jt=p, the transition proba¬ 
bilities to tilde and non-tilde phonon states are connected 
as 


$ 


(T) 

Ji 



( 10 ) 


In [15], a similar relation is obtained for charge-neutral 
transitions and it is referred to as the principle of de¬ 
tailed balance. However, in contrast to the case of charge- 
neutral transitions, the detailed balance relation (10) in¬ 
cludes a phonon energy rather than a transition energy. 
The latter is the energy transferred to the parent nucleus 
and for charge-changing transitions is given by [25] 


Eji — JJ Ji ~F &np • 

Elf = -4f. (11) 


Here and below, the upper (lower) sign refers to the neu¬ 
trino (antineutrino) cross section. In the above expres¬ 
sion, Gf is the Fermi constant for the weak interaction, 
9c is the Cabbibo angle, and E e and p e are the energy 
and momentum of the outgoing electron or positron. 

In Eq. (13), the contributions, <Tq L , for the Coulomb 
and longitudinal components, and for the transverse 
electric and magnetic components, are written as 


a CL = (1 + acosO)|(Ji||Afj||0(T))| 2 + (1 + cos© 


— 25sin 2 0)|(Ji||Lj||O(T))| s 


Ejj 

L ~ 


(1 + COS0) + c 


x 2Re(Ji\\Lj\\0(T))(Ji\\Mj\\0(T))*, (14) 


cj j t = (1 — acos© + 6 sin 2 0) |j(Ji||T“ ag ||0(T , ))| 2 + 
r E v + e , 


KJi||7?||0(T)>|‘ 


T 


-(1 — acos0) — c 


x 2 Re( Ji||Tj nag ||0(T))( Ji||Tj ||0(T))*, (15) 


where S np = AX np + A M np , and AX np is the difference 
between neutron and proton chemical potentials in the 
parent nucleus and A M np = 1.293 MeV is the neutron- 
proton mass splitting. It is obvious that non-tilde (tilde) 
charge-changing phonon states do not necessarily corre¬ 
spond to upward (downward) transitions: due to S np , 
some non-tilde (tilde) phonon states may have negative 
(positive) transition energy. With Eq. (11), relation (10) 
is rewritten as 

= exp(- ^ ) r T ^ )^t ) - ( 12 ) 

Thus, for each n —> p (p —> n) upward transition with 
energy E there is an inverse downward transition p —> n 
(n —>• p) with energy —E and the respective transition 
probabilities are connected by (12). In the Appendix we 
show that in this form the detailed balance for charge¬ 
changing transitions in the ensemble of hot nuclei can be 
derived in a model independent way. 

B. Cross sections in supernova environments 

In the derivation of the (anti)neutrino absorption cross 
section for a hot nucleus under supernova conditions we 
follow the Walecka-Donnelly formalism [26, 27] which is 
based on the standard current-current form of the weak 
interaction Hamiltonian. After a multipole expansion of 
the weak nuclear current the temperature-dependent dif¬ 
ferential cross section for a transition from the thermal 


where 0 is the lepton scattering angle and the nota¬ 
tion |Ji) is used to denote both the non-tilde and the 
tilde thermal one-phonon states. In the former case the 
transition energy Eji = Ejf\ while in the latter case 
Eji = Ejp. Thus, the energy of the outgoing lepton is 
E l e = Ev — Eji. The parameters a, b and c are obtained 
from the relations 



and the absolute value of the three-momentum transfer 
q is given 

q = ^E 2 Ji +2E e E„{l-acosQ) - (?n e c 2 ) 2 . (17) 

The multipole operators Mj , Lj, jf, and J“ ag denote 
the charge, longitudinal, and transverse electric and mag¬ 
netic parts of the hadronic current, respectively, as de- 
finedvin [26, 27]. For the vector, axial-vector, and pseu¬ 
doscalar form-factors which describebthe internal struc¬ 
ture of the nucleon we use parametrization from Ref. [28] 
(see also Ref. [29]). 

For charged-current reactions, the cross section (13) 
must be corrected for the distortion of the outgoing lep¬ 
ton wave function by the Coulomb field of the residual 
(daughter) nucleus. The cross section can either be mul¬ 
tiplied a Fermi function F(Z,E) (see, e.g., Ref. [30]), 
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or, at higher energies, the effect of the Coulomb field 
can be described by the effective momentum approxima¬ 
tion (EMA) [31]. In the present study, by following the 
prescription from [32, 33], we choose an energy point in 
which both approaches predict the same values. Then 
the Fermi function is used below the point and the EMA 
is adopted above it. 

Furthermore, in the supernova environment nuclei are 
surrounded by a nearly degenerate electron gas. Thus, 
neutrino absorption can be strongly reduced by Pauli 
blocking of the final electron states. The blocking factor 
(1 — f(E e )) in Eq. (13) accounts for this effects, where 
f(E e ) is the Fermi-Dirac distribution with temperature 
T and the chemical potential /i e - . The positron distribu¬ 
tion is defined in the same way with fi e + = — fi e - . There¬ 
fore, we can neglect the blocking factor for antineutrino 
absorption. 

The total cross section a(E V ,T), as a function of tem¬ 
perature and incoming (anti)neutrino energy, is obtained 
from the differential cross sections (13) by integrating 
over the scattering angle and summing over all possible 
final thermal one-phonon states 

<j(E„,T) = 2-tt d cos 0 

=a en (E v ,T) + a ex (E v ,T). (18) 

Here, the total cross section is split into two parts: 
a en (E u ,T) describes the endoergic neutrino absorption 
and includes only upward transitions (Eji > 0), while 
a ex (E u , T) corresponds to the exoergic process associated 
with downward transitions [Eji < 0). The contribution 
a ex (E v , T) dominates the cross section for vanishing neu¬ 
trino energies, E„ k, 0. It is apparent that for beta-stable 
nuclei the exoergic absorption is only possible at finite 
temperatures and due to transitions from thermally ex¬ 
cited states. 

For low-energy (anti)neutrinos, i.e., in the long wave¬ 
length limit q —» 0, the allowed 1 + multipole operator 
reduces to the Gamow-Teller form 


GT T =g A at T , (19) 

where gA = —1.2599 [34] is the axial weak coupling con¬ 
stant. Considering only Gamow-Teller transitions and 
taking into account detailed balance according to (10), 
the cross section (18) can be written as 


<t{E v ,T) = 


(Gp cos 9c ) 2 


x - E^ ) ) 2 ^ ) F(±Z + 1, E e ) 

i 

+ £(£„ + E^f exp(-^)s *?>F{±Z + 1, E e )} 


( 20 ) 


Here for simplicity the electron rest mass has been ne¬ 
glected, i.e. E e ss p e , and the blocking factor for the 


outgoing lepton has been dropped. The matrix elements 
denote the transition strength of the GT^ opera¬ 
tor. The first term in Eq. (20) implies summation over 
non-tilde 1 + states. Within the shell model approach 
this term reduces to the ground-state contribution (see 
Eq. (2) in Ref. [7]). However, in the TQRPA this term 
appears to be temperature dependent due to violation of 
Brink’s hypothesis. The second term accounts for tran¬ 
sitions to tilde states and it is an analog of the back- 
resonance contribution within the shell-model approach. 
Here, we would like to stress that we cannot associate 
transitions to non-tilde (tilde) states with the endoergic 
(exoergic) absorption. As was mentioned above, due to 
S np some non-tilde (tilde) states may correspond to down¬ 
ward (upward) transitions and, hence, contribute to a ex 
(u e n ) ■ 


III. RESULTS AND DISCUSSION 


We employ the theoretical framework presented above 
to study thermal effects on the (v e ,e~) and (77 e , e + ) cross 
sections for two sample nuclei, 56 Fe and 82 Ge. The iron 
isotope is among the most abundant nuclei at the early 
stage of the core-collapse, while the neutron-rich germa¬ 
nium isotope can be considered as the average nucleus at 
later stages [35]. 

To describe charge-changing transitions in a hot nu¬ 
cleus we use the same phenomenological nuclear Hamil¬ 
tonian as in [18]. The Hamiltonian consists of spherically 
symmetric Woods-Saxon potentials for protons and neu¬ 
trons, BCS pairing interactions, and separable multipole 
and spin-nrultipole residuals interaction in the particle- 
hole channel. We neglect particle-particle interactions 
except for the BCS pairing forces. This Hamiltonian 
is usually referred to as the quasiparticle-phonon model 
(QPM) [36]. For the two nuclei considered, the param¬ 
eters of the QPM Hamiltonian are fixed in the same 
manner as in [15, 18]. Here we just mention that solv¬ 
ing the BCS equations at zero temperature we get the 
following proton and neutron pairing gaps: A p ( n ) = 
1.57(1.36) MeV for 56 Fe and A p(n) = 1.22(0.0) MeV for 
82 Ge. Thus, the critical temperature (T cr « 0.5A r ) when 
the pairing phase transition occurs is T cr « 0.8 MeV for 
the iron isotope and T CI « 0.6 MeV for the germanium 
isotope. 

A separable form of the residual interaction allows us 
to reduce the TQRPA equations to a secular equation. 
The explicit form for charge-changing excitations as well 
as expressions for transition strengths (9) are given in [17, 
18]. Moreover, it has been proved in Ref. [17] that the 
TQRPA fulfills the Ikeda sum rule for allowed Fermi and 
GT transitions. 
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FIG. 1. (Color online) Cross sections for neutrino (upper panels) and antineutrino (lower panels) absorption reactions for 56 Fe 
and 82 Ge. The total cross sections (solid lines) include contributions of the J n = 0 ± — S 1 * 1 multipoles. The dashed lines show 
the cross sections without the 1 + contribution. The dash-dotted lines correspond to the 1 + contributions calculated with the 
full momentum-dependent transition operator whereas the 1 + contributions calculated with the Gamow-Teller operator (19) 
are shown by the dash-double dotted lines. 


A. Ground-state cross sections 

As in [15], before proceeding to discuss thermal ef¬ 
fects on neutrino-nucleus absorption cross sections, we 
consider the ground-state (T = 0) cross sections for 
56 Fe(i/ e , e _ ), 56 F e(P e , e + ), 82 Ge(^ e , e ~), and 82 Ge(F e , e + ) 
reactions and compare the our results with those avail¬ 
able from other theoretical studies. In Fig. 1, the cal¬ 
culated total cross sections are shown as functions of 
the incident (anti)neutrino energy E v . We find that for 
E v < 100 MeV the reactions considered are dominated 
by the multipole transitions J 7r < 3^, while contribu¬ 
tions from higher multipoles are only a few percent of 
the total cross sections. 

For 56 Fe the cross sections increase sharply as E u ap¬ 
proaches the reaction threshold Q_ = 4.56 MeV for 
j/ e -absorption and Q+ = 4.71 MeV for 77 e -absorption 
(Q^ = Mf + m e c 2 — Mi, where Mi j are the masses of the 
parent and daughter nuclei). As the neutron number in¬ 
creases, the threshold energy for i/ e -absorption decreases. 
For the neutron-rich nucleus 82 Ge the Q_-value becomes 
negative (Q_ = —4.71 MeV) allowing neutrino absorp¬ 
tion for all energies. Contrary to this, for the (F e ,e + ) 
reactions the Q + -values become less favorable with in¬ 
creasing neutron excess. The antineutrino has to over¬ 
come a noticeable threshold (Q + = 13.58 MeV) to be 
absorbed by 82 Ge and the corresponding cross section is 
considerably lower in comparison to the neutrino one. 


In addition to the total cross sections, contributions 
from the 1 + multipole channel are shown in Fig. 1. 
Referring to the plots in Fig. 1 for the 56 F e{y e ,e~), 
56 Fe(V e , e + ), 82 Ge(r' e , e~) reactions the cross sections are 
dominated by allowed 1 + transitions for energies up to 
E„ = 30 MeV, with contributions from other multipoles 
being much smaller. In the 82 Ge(y e ,e + ) cross section, 
however, the 1 + contribution is negligible and forbidden 
transitions dominate the reaction. We find that 1“ and 
2“ transitions mainly contribute to the cross section for 
E v < 50 MeV. 

In Fig. 1 we also analyze the effect of the full in¬ 
dependent 1 + transition operator instead of its long- 
wavelength limit. As was mentioned above, in the latter 
case the 1 + transition operator reduces to the Gamow- 
Teller operator (19) and the ground-state cross section 
is given by the first term in Eq. (20). It should be 
noted that to make a comparison with the shell-model 
calculations [37, 38] more transparent, we use the same 
quenching factor for the axial weak coupling constant, 
g* A = 0.74^. By comparing in Fig. 1 the 1 + and Gamow- 
Teller contributions to the cross sections, we conclude 
that for energies E„ < 30 MeV the application of the 
GT operators instead of the g-dependent 1 + operator 
is fully justified. Therefore, we conclude that the low- 
energy ground-state cross sections for the 56 F e(u e ,e~), 
56 Fe(u e , e + ), 82 Ge(iz e ,e _ ) reactions are completely dom¬ 
inated by GT transitions. In contrast, all GT+ transi- 
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TABLE I. Total 56 F 'e(v e ,e~) cross sections for selected neu¬ 
trino energies E„. The present QRPA results (second column) 
are compared with those from [39] and with the hybrid ap¬ 
proach results [37]. The cross sections are given in units of 
1CP 42 cm 2 , with exponents given in parentheses. 


E v (MeV) 

QRPA 

SkQRPA [39] 

Hybrid[37] 

10 

2.39(—1) 

3.63(+0) 

6.61(—1) 

15 

2.35(+0) 

1.73(+1) 

6.45(+0) 

20 

1.91(+1) 

5.26(+l) 

2.93(+l) 

25 

6.19(+1) 

1.25(+2) 

7.33(+l) 

30 

1.34(+2) 

2.33(+2) 

1.40(+2) 

40 

3.88(+2) 

5.44(+2) 

3.71(+2) 

50 

8.47(+2) 

9.83(+2) 

7.98(+2) 

60 

1.58(+3) 

1.67(+3) 

1.38(+3) 

70 

2.61(+3) 

2.59(+3) 

2.42(+3) 

80 

3.95(+3) 

3.73(+3) 

3.60(+3) 

90 

5.53(+3) 

5.07(+3) 

4.98(+3) 

100 

7.26(+3) 

6.60(+3) 

6.52(+3) 


TABLE II. Same as Table I, but for the 56 Fe(^ e , e + ) reaction. 


E v (MeV) 

QRPA 

SkQRPA [39] 

Hybrid [38] 

10.0 

4.92(—1) 

2.95(+0) 

0(+0) 

12.5 

1.34 (+0) 

6.09(+0) 

0(+0) 

15.0 

2.59(+0) 

1.03(+1) 

2(+0) 

20.0 

6.48(+0) 

2.17(+1) 

7(+0) 

25.0 

1.28(+1) 

3.74(+l) 

1.6(+1) 

30.0 

2.26(+l) 

5.74(+l) 

3.0(+l) 

40.0 

5.71(+1) 

1.20(+2) 

8.4(+l) 

50.0 

1.17(+2) 

2.09(+2) 

1.81(+2) 


tions of are essentially blocked in 82 Ge. For nuclei with 
Z < 40 and IV > 40 blocking occurs because the valence 
protons are in the pf shell, while the valence neutrons oc¬ 
cupy already the next major shell (sdg shell). In the next 
section it will be shown, however, that thermal effects un¬ 
block GT + transitions in 82 Ge and for typical supernova 
temperatures the low-energy 82 Ge(F e ,e + ) cross sections 
are also dominated by the GT contributions. 

In Tables I and II, we compare the calculated ground- 
state cross sections for (anti)neutrino absorption by 56 Fe 
to those obtained from self-consistent QRPA calcula¬ 
tions with Skyrrne forces [39] and with the hybrid ap¬ 
proach (large-scale shell-model calculations for GT con¬ 
tributions and RPA for other multipoles) [37, 38]. As 
is seen, the three models, although based on different 
microscopic pictures, predict rather similar neutrino ab¬ 
sorption cross sections. We note, however, that for low 
energies (E„ < 40 MeV) the present results are closer to 
the hybrid approach results than those of Ref. [39]. It is 
evident that the reason why the QRPA and the hybrid 
approach cross sections are systematically lower than 
those of the Skyrrne based calculations at E v < 40 MeV 



FIG. 2. Contribution of the multipole transitions J n = 
(V — 3* to the cross section for the S6 Fe(^ e , e“) 56 Co reac¬ 
tion at E v = 40, 60, 80 MeV. The present QRPA results 
are compared with those obtained from RQRPA [41] and 
SkQRPA [39] calculations. 


is caused by differences in the GT_ strength distribu¬ 
tions. Although the authors of Ref. [39] do not provide 
the GT_ distribution in 56 Fe, it seems that the Skyrrne 
based calculations result in a larger strength below the 
GT_ resonance compared to the other two approaches. 
Other possible reasons for the discrepancy could be a 
lower energy of the GT_ resonance and a larger total GT 
strength obtained from the calculations with the Skyrrne 
interaction. 

Experimental results for neutrino absorption cross sec¬ 
tions are rather limited. The KARMEN Collaboration 
has measured the flux averaged 56 Fe(r' e , e _ ) 56 Co cross 
section for the neutrino spectrum from the muon decay 
at rest and obtains (a) = [251 ± 83(stat.) ± 42(syst.)] x 
10 -42 cm 2 [40]. Our result, (a) = 223 x 10 _42 cm 2 , 
is in excellent agreement with the experimental value. 
Note, that the hybrid approach predicts (a) = 240 x 
10 _42 cm 2 [37], while the QRPA calculations with Skyrrne 
forces gives (a) = 352 x 10 _42 cm 2 [39]. 

For 56 Fe a detailed demonstration of the most impor¬ 
tant multipole contributions to neutrino absorption cross 
section is presented in Fig. 2 at three neutrino energies, 
E u = 40, 60, and 80 MeV. As expected, at relatively low 
neutrino energies ( E „ < 40 MeV) the dominant contribu¬ 
tion to the cross section originates from 1 + transitions. 
With increasing E v , however, contributions from other 
multipole transitions become important. In particular, 
at E„ = 80 MeV, the dominant contribution comes from 
1“ transitions, but other multipoles, e.g., 1 + , 2 _ , 2 + , 
and 3 + , also play an important role. 

In Fig. 2, we also compare the calculated multipole 




decomposition of the cross section for 56 Fe with those 
from the SkQRPA calculations [39] and the relativistic 
QRPA (RQRPA) calculations [41]. As is evident from the 
figure, the latter model predicts somewhat smaller cross 
section at high neutrino energies (E v > 60 MeV), whereas 
we observe an excellent agreement between the results 
of the present QRPA and the SkQRPA. Specifically, in 
accordance with Ref. [39, 41], we find that 0 + allowed 
transitions only marginally contribute to the 56 Fe(i/ e , e~) 
reaction. This finding is true for other three reactions 
and for finite temperatures as well. For this reason, in 
the discussion below, we will always imply 1 + multipole 
channels when referring the allowed transitions. 


B. Finite temperatures 

In discussing v e and v e absorption reactions under 
supernova conditions we will follow the line of our re¬ 
cent work [15] and consider first thermal effects on the 
strength distribution of GT transitions which dominate 
the reactions at low energies. 

In Fig. 3, we display on a logarithmic scale the GT_ 
and GT + distributions in 56 Fe calculated for the ground- 
state (T = 0) and at three stellar temperatures above 
the critical one, T = 0.86 MeV (10 10 K), 1.29 MeV (1.5 x 
10 lo K), and 1.72MeV (2 x 10 lo K). These temperatures 
roughly correspond to three stages in the collapse evolu¬ 
tion [9]. We emphasize that the distributions are plot¬ 
ted with the bare GT operators at ^ as functions of the 
transition energy (11). For each temperature we show 
the value of S np in Eq. (11) as well as the total transi¬ 
tion strengths SL and S+. Note that S- and S+ values 
calculated for the ground-state satisfy the Ikeda sum rule 
S- — S+ = 3 (N—Z) (a small deviation is caused by the in¬ 
completeness of our single-particle basis) but noticeably 
overestimate the experimental data ( S _ = 9.9 ± 2.4 [42], 
S+ = 2.9 ± 0.3 [43]). This overestimation is common for 
any QRPA calculations of GT strength and is remedied 
by a quenched value for the axial weak coupling con¬ 
stant gA■ One might notice that our QRPA calculations 
fairly well reproduce the experimental centroid energies 
for both GT+ [43] and GT_ [42] distributions in 56 Fe. In 
this respect the present calculations are consistent with 
the large-scale shell-model calculations [44]. Of course, 
in contrast to the LSSM approach, the QRPA cannot 
recover all nuclear correlations needed to reproduce the 
fragmentation of the GT strength. 

For 56 Fe, due to a relatively small absolute value 
of S np , the major part of the upward (downward) 
GT strength corresponds to transitions to non-tilde 
(tilde) phonon states. According to our QRPA calcu¬ 
lations, the principal contribution to the GT_ (GT + ) 
strength in 56 Fe comes from the I/7/2 —> I/5/2 neutron- 
to-proton (proton-to-neutron) single-particle transition 
which forms a resonance peak at E « 14 MeV (« 6 MeV). 
Because the Brink hypothesis is not valid within the 
TQRPA, we observe an evolution of the GT resonances 



FIG. 3. (Color online) Temperature evolution of GT_ (left 
panels) and GT+ (right panels) strength distributions for 56 Fe 
vs transition energy. The solid (dashed) lines refer to tran¬ 
sitions to non-tilde (tilde) thermal one-phonon states. The 
arrows indicate the ground-state reaction thresholds for neu¬ 
trino ( Q- = 4.56 MeV) and antineutrino (Q+ = 4.71 MeV) 
absorption. 


with temperature. Namely, when the temperature is in¬ 
creased to 0.86 MeV, the transition energy is lowered by 
~ 1.2 MeV for the GT_ resonance and ~ 2.1 MeV for the 
GT + resonance. This decrease in energy is mainly at¬ 
tributed to the vanishing of pairing correlations: at tem¬ 
peratures above the critical one no extra energy is needed 
to break a proton (neutron) Cooper pair and as a conse¬ 
quence the GT + (GT_) resonance moves to lower ener¬ 
gies. However, not only pairing effects lead to downward 
shift of the GT resonances. It can be easily seen from 
the structure of the TQRPA equations [17] that quasi¬ 
particle thermal occupations factors will appear which 
screen the interaction term. Due to the thermal block¬ 
ing of the proton-neutron repulsive residual interaction, 
a further increase in temperature to 1.72 MeV decreases 
the GT_ and GT+ resonances in 56 Fe by ~ 0.3 MeV 
and ~ 0.5 MeV, respectively. As was mentioned in the 
introduction, the observed downward shift of the GT 
strength is not present in large-scale shell-model calcula¬ 
tions which are partially based on Brink’s hypothesis. In 
contrast, the finite-temperature relativistic QRPA calcu¬ 
lations in Ref. [45] and shell-model Monte-Carlo calcula¬ 
tions in Ref. [10] show similar features for the changes of 
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the GT resonance energy. 

Finite temperature also affects the low-energy GT_ 
strength in 56 Fe: due to the gradual disappearance of 
pairing it partially shifts below the ground-state reac¬ 
tion threshold and becomes more fragmented. The frag¬ 
mentation arises from the thermal smearing of the nu¬ 
clear Fermi surface that unblocks low-energy particle- 
particle and hole-hole GT_ transitions. Here, particle 
(hole) denotes a state above (below) the Fermi level. 
Moreover, in accordance with detailed balance (12), the 
temperature rise exponentially increases the strength of 
negative-energy downward transitions. Referring to the 
figure, the dominant part of the downward GT_ strength 
around E ss —4.1 MeV originates from a transition in¬ 
verse to the GT + resonance, that is due to the transi¬ 
tion from the thermally populated neutron orbit I/5/2 
to the lower lying proton orbit 1 / 7 / 2 - In contrast, tran¬ 
sitions inverse to the GT_ resonance are suppressed by 
the energy-dependent exponent in Eq. (12) and the GT + 
downward strength is dominated by transitions inverse 
to low-energy GT_ transitions. It should be mentioned 
that although thermal effects unblock some new GT tran¬ 
sitions, the total GTip strength in 56 Fe slightly decreases 
with temperature. Nevertheless, the TQRPA preserves 
the Ikeda sum rule. 

It is obvious that the violation of Brink’s hypothesis 
should affect the downward GT strength. In [15], this 
influence was studied thoroughly for charge-neutral GT 
transitions in 56 Fe by comparing the GT running sums 
derived with and without Brink’s hypothesis. Applying 
the same method, we find that both the shift of the GT+ 
resonance to lower energies and the thermal unblocking 
of low-lying GT_ strength enhance the GT_ and GT+ 
downward strengths in 56 Fe. 

Fig. 4 shows the temperature evolution of the GT 
strength function for the neutron-rich nucleus 82 Ge. The 
distributions are displayed at the same temperatures as 
those for 56 Fe in Fig. 3. As evident from the figure, some 
GT strength associated with non-tilde (tilde) phonon 
states appears at negative (positive) transition energies 
which is due to a large value of S np . The temperature 
rise affects the GT_ and GT + distributions in a different 
way. For the GT_ distribution, a temperature increase 
essentially affects only the downward part, yielding some 
new strength below —10 MeV. By comparing the left and 
right panels in Fig. 4, we conclude that this thermally 
unblocked strength corresponds to transitions inverse to 
the GT + resonance. The GT_ resonance is weakly sen¬ 
sitive to thermal effects. By increasing the temperature 
to T = 1.72 MeV, the excitation energy gets only low¬ 
ered by ~ 0.5 MeV. Since there are no neutron pairing 
correlations in 82 Ge, this lowering is solely caused by the 
softening of the residual interaction. 

In contrast, thermal effect are significant for the GT+ 
resonance and they clearly demonstrate a violation of the 
Brink hypothesis within the TQRPA. Namely, with ris¬ 
ing temperature the resonance peak shifts by ~ 7 MeV to 
lower energies and the total GT + strength shows a non- 



E (MeV) 


FIG. 4. (Color online) Same as Fig. 3 but for 82 Ge. Due to 
Pauli blocking of allowed p —>• n transitions, the total GT_ 
strength is many times larger than the GT+ strength and 
S- « 3 (N -Z) = 54. 

monotomic dependence: after an initial decrease, it grad¬ 
ually increases. In [18], a detailed investigation within 
the TQRPA approach was performed of thermal effects 
on the GT + strength in the neutron-rich nucleus 76 Ge. 
The conclusions remain valid for 82 Ge, as well. Briefly, 
for neutron-rich nuclei with TV > 40 and Z < 40 the Inde¬ 
pendent Particle Model predicts that all GT+ transitions 
of valence protons are Pauli blocked due to the complete 
occupation of the pf neutron orbits. However, both ther¬ 
mal excitations and pairing correlations promote protons 
to the sdg shell and make possible lffg/ 2 —> lg", 2 particle- 

particle transition in 82 Ge possible. The important point 
is that the transition energy drastically decreases with 
temperature. Namely, at T < T cr , when pairing cor¬ 
relations are important, the transition energy is given 
by e lg n /2 + s lg p /2 + S pn » 20 MeV, while at T > T cr , 
when thermal effects become significant, it is given by 
£is" /2 ~ £1 g p /2 + Spn « 13 MeV. In addition, the transi¬ 
tion strength below the GT + resonance becomes increas¬ 
ingly unblocked with temperature. In particular, a strong 
peak due to the l/^ 2 —> 1 /"^ 2 hole-hole transition ap¬ 
pears at E ss 11 MeV. Both unblocking mechanisms are 
severely suppressed in the vicinity of the critical tem¬ 
perature ( T CI ft: 0.6 MeV in 82 Ge), i.e. when pairing 
correlations vanish while thermal effects are not yet suf- 
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FIG. 5. (Color online) (Upper panels) Total neutrino ab¬ 
sorption cross sections for 56 Fe and 82 Ge at three different 
temperatures relevant for core collapse. For comparison, the 
ground-state cross sections are also shown. (Middle panels) 
Contributions from forbidden transitions to the cross sections. 
(Lower panels) Temperature dependence of the ratio of exo- 
ergic absorption to the reaction cross section. 


ficiently strong to significantly occupy the lg g / 2 proton 
orbit or unblock the I/ 5/2 neutron orbit. Therefore, the 
total GT_|_ strength decreases at T w T cr . The ther¬ 
mal effects on the GT + strength in neutron-rich nuclei 
discussed here were predicted in [35] and also confirmed 
in [45] by calculations based on the finite-temperature 
relativistic QRPA. We also note that, due to the large 
S np , downward GT + transitions appears to be highly sup¬ 
pressed in 82 Ge. 

Let us now illustrate the influence of the thermal ef¬ 
fects discussed above on the neutrino and antineutrino 
absorption by 56 Fe and 82 Ge. In Fig. 5, we compare 
the (u e ,e~) ground-state cross sections with those cal¬ 
culated at the three core-collapse temperatures. For en¬ 
ergies E v < 30 MeV, which are typical for supernova 
neutrinos, the cross sections are dominated by allowed 
GT_ transitions at all temperatures. This is verified in 
the middle panels of Fig. 5, where the overall contribu¬ 
tion of forbidden transitions to the finite temperature 
cross sections is given. Although the forbidden cross sec¬ 
tion increases with temperature, it contributes less than 
30% all the way up to E v = 30 MeV. The lower pan¬ 
els of Fig. 5 show the ratio of exoergic absorption to the 
reaction cross section 


P(E*,T) 


a ex (E„,T) 
<t{E v ,T) ’ 


( 21 ) 


where a ex (E v ,T) only accounts for negative-energy 


downward transitions induced by neutrino absorption. 

When considering the 56 Fe(i/ e , e _ ) cross section, we 
observe that thermal effects are unimportant for E v > 
20 MeV, i.e. when incoming neutrinos have sufficiently 
large energy to excite the GT_ resonance. Note, that 
a downward shift of the GT_ resonance only marginally 
affects the cross section at such high neutrino energies. 
Thermal effects become pronounced though for lower 
neutrino energies since finite temperature removes the 
reaction threshold and significantly enhances the cross 
section. For energies E v < 5 MeV, the ratio /3 > 0.5 and, 
hence, the observed enhancement is mostly caused by 
downward GT_ transitions from nuclear excited states. 
Moreover, for E v « 0 the thermally unblocked downward 
transitions completely dominate the reaction, increasing 
the cross section by more than an order of magnitude 
when the temperature rises from 0.86 MeV to 1.72 MeV. 
The role of the exoergic absorption decreases with in¬ 
creasing neutrino energy and for intermediate energies, 
5 MeV < E v < 20 MeV, thermal effects on the GT_ res¬ 
onance and its low-energy tail become important for the 
cross section enhancement. 

Referring to the right panels in Fig. 5 it is shown that 
the 82 Ge(i/ e ,e _ ) cross section to a much lesser extent 
depends on temperature than that for the 56 Fe(^ e ,e“) 
reaction. This result can be understood as follows. In 
82 Ge, due to the negative Q_-value, the downward GT_ 
transitions dominate the ground-state reaction up to 
E v ss 15 MeV. As discussed above, the GT_ distribution 
in 82 Ge is little affected by the temperature rise, which 
merely causes some additional strength at E < —10 MeV. 
This thermally unblocked downward strength becomes 
competitive with the ground-state contribution only at 
rather low energies (E u < 5 MeV). As a result, thermal 
effects on the 82 Ge(^ e , e ~) reaction are noticeably milder 
than in the previous case and a temperature rise from 
T = 0 to 1.72 MeV enhances the low-energy cross section 
only by a factor of about four. This observation is in 
line with Ref. [7], where it was shown that the thermal 
enhancement of the neutrino absorption cross section is 
reduced the more neutron-rich the nucleus is. 

The results for 77 e -absorption are shown in Fig. 6 . For 
the 56 Fe(z7 e , e + ) reaction the cross section demonstrates 
the same trend as discussed above for the neutrino ab¬ 
sorption by 56 Fe: (i) The gap in the cross section dis¬ 
appears and the low-energy cross section increases with 
temperature, (ii) with increasing E„ the cross sections 
at different temperatures converge, i.e., thermal effects 
become unimportant. From the lower-left panel we con¬ 
clude that a significant enhancement of the low-energy 
cross section relative to ground-state calculations comes 
from the increasing contribution of downward GT_ tran¬ 
sitions from nuclear excited states. The ratio j3 gradually 
reduces at E v > 5 MeV, since at those energies the exci¬ 
tation of the GT + resonance becomes possible. However, 
at T = 1.72 MeV the exoergic component of the cross sec¬ 
tion appears to be comparable with the endoergic one up 
to E v « 15 MeV. 
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FIG. 6. (Color online) The same as in Fig. 5, but for the FIG. 8. (Color online) The same as in Fig. 5, but now includ- 
antineutrino absorption reaction. ing the outgoing electron blocking 
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FIG. 7. (Color online) Relative contribution ai/a of the ith 
negative-energy GT state to the finite-temperature cross sec¬ 
tion at E v = 0. 


It should be noted that although the calculated 
56 Fe(i7 e ,e“) and 56 Fe(77 e ,e + ) cross sections show the 
same trend, the observed thermal enhancement at low 
energies is caused by different types of downward tran¬ 
sitions. To show that, we have calculated the relative 
contribution eq/c t of different negative-energy GT states 
to the cross section at E v = 0. The results are depicted 
in Fig. 7. As indicated in the figure, the i/-absorption is 
completely dominated by the thermally unblocked transi¬ 
tion inverse to the GT + resonance. For the 77-absorption, 
however, the main contribution to the reaction is given 
by transitions inverse to the GT_ low-energy strength, 
while the "inverse" GT_ resonance is suppressed by the 


Boltzmann factor in Eq. (20). Only at T = 1.72 MeV, the 
"inverse" GT_ resonance gives a noticeable contribution 
to the absorption of low-energy antineutrinos. 

As illustrated in the right panels of Fig. 6, the low- 
energy 82 Ge(77 e ,e + ) cross section at finite temperature 
is also predominantly mediated by GT+ transitions. Al¬ 
though the cross section remains substantially smaller as 
compared with the other three reactions considered, it 
demonstrates a strong temperature dependence for an¬ 
tineutrinos with Ev < 15 MeV, i.e. below the ground- 
state reaction threshold. Since the downward GT + tran¬ 
sitions are suppressed in 82 Ge, the cross section en¬ 
hancement reflects the thermal unblocking of the up¬ 
ward GT + strength. Because of the thermal unblocking, 
the energy below which the GT + contribution is more 
than half of the total cross section shifts to higher val¬ 
ues: at T = 0.86 MeV this energy about 15 MeV, at 
T = 1.29 MeV it is about 20 MeV, and at T = 1.72 MeV 
it is about 25 MeV. For higher energies forbidden transi¬ 
tions dominate the process like at T = 0 and the cross- 
section depends only weakly on temperature. 

As was mentioned above, during the core contraction 
neutrino absorption by nuclei is hindered by Pauli block¬ 
ing of the phase space for the outgoing electron. To study 
this effect within the TQRPA, following Ref. [7], we have 
calculated the ^-absorption cross sections by introduc¬ 
ing a blocking factor (1 — f(E e )) defined at three dif¬ 
ferent sets of temperature and chemical potential (/i e in 
MeV): (T,p e ) = (0.86, 8.3), (1.29,18.1), and (1.72,36.2). 
The results are shown in Fig. 8 for the 56 Fe(i/ e ,e _ ) and 
82 Ge(i/ e ,e _ ) reactions. As can be seen, the neutrino ab¬ 
sorption cross sections are drastically reduced due to elec¬ 
tron blocking in the final state. Moreover, as the chem- 
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FIG. 9. (Color online) Comparison of the cross sections of 
neutrino absorption on the hot nucleus 56 Fe calculated within 
the TQRPA and the LSSM approach (see Ref. [7], Fig. 1). 
The cross sections in the upper (lower) panel are derived with¬ 
out (with) the final-state electron blocking. The respective 
temperatures and chemical potentials are given in MeV. 


ical potential increases faster than the temperature the 
cross sections decrease with temperature. It is obvious 
that the absorbtion due to de-excitation of thermally ex¬ 
cited states is less affected by the Pauli blocking since the 
outgoing electron gains energy. This is clearly demon¬ 
strated by the plots in the lower panels of Fig. 8 where 
the relative contribution of negative-energy transitions is 
shown. We thus conclude that owing to blocking effect 
negative energy transitions dominate the cross sections 
up to neutrino energies E v > /j, e . Another consequence of 
the blocking is that downward 1 hui forbidden transitions 
become important with increasing /i e . Without Pauli 
blocking their contribution in negligible due to a small 
Boltzmann weight. However, as shown in the middle 
panels of Fig. 8, for (T,/z e ) = (1.29,18.1) and (1.72,36.2) 
their contribution is comparable or even larger than those 
of the allowed transitions. According to our calculations 
this downward forbidden contribution is mainly due to 
CT, 1 _ ,2 _ transitions. 

In Fig. 9, the 56 Fe(i/ e ,e _ ) cross sections are displayed 
together with the results of the large-scale shell-model 
calculations (see Fig. 1 in Ref. [7]). Even if the general 
behaviour of the cross sections as functions of the neu¬ 


trino energy and temperature is in agreement in both ap¬ 
proaches, the TQRPA results are noticeable larger than 
the LSSM ones but the discrepancy reduces with increas¬ 
ing neutrino energies. To understand the cause of the 
discrepancy, let us consider first the cross sections cal¬ 
culated neglecting the final-state electron blocking. For 
E u < 5 MeV the TQRPA cross sections exceed the LSSM 
values by two to three orders of magnitude. For such 
low energies, in both approaches, the neutrino absorp¬ 
tion on 56 Fe is dominated by GT_ downward transi¬ 
tions from thermally excited states. The differences in 
the description of such transitions explain the differences 
between the TQRPA and shell-model results. As the 
TQRPA is based on the grand canonical ensemble, the 
upward and downward GT± strength are connected by 
the detailed balance. As shown in Fig. 7, the absorp¬ 
tion of low-energy neutrinos on 56 Fe is dominated by the 
downward GT_ transition inverse to the GT + resonance. 
This corresponds to an excitation energy u>i ~ 6 MeV 
in the Boltzmann factor of Eq. (20). Within the shell- 
model calculations, downward transitions are included by 
back-resonances, that is by inverting the GT + strength 
distribution of the daughter nucleus. In [7], the back- 
resonances are built on the lowest states of 56 Co and the 
bulk of the back-resonance strength in 56 Fe is located 
at an excitation energy of Ei « 7 — 9 MeV. This ex¬ 
citation energy is somewhat higher than Wj. Moreover, 
within the TQRPA the downward strength concentrates 
in a single state, while the shell model GT + strength 
for 56 Co is highly fragmented owing to multinucleon cor¬ 
relations [46]. It is clear that both these factors sup¬ 
press the contribution of downward transitions within the 
LSSM. To see whether the TQRPA reliably predicts the 
strength of negative-energy transitions, one must go be¬ 
yond the TQRPA. For a separable residual interaction 
used here this can be done following the method devel¬ 
oped within the QPM, that is by taking phonon cou¬ 
pling into account. On the other hand we should note 
that due to violation of Brink’s hypothesis some back- 
resonances built on high-lying excited states of 56 Co may 
be located at the same energies as those built on the nu¬ 
clear ground and low-lying states. Due to an increasing 
density of states, the contribution of such back-resonance 
states may be substantial and, therefore, their inclusion 
into the shell-model calculations may improve the agree¬ 
ment between the TQRPA and shell-model results. For 
5 < £7^ <15 MeV, thermally unblocked low-energy GT_ 
transitions come into play, and since such transitions do 
not appear within the shell-model based calculation they 
also cause the excess of the TQRPA cross sections over 
the LSSM ones. At E v > 15 MeV, jz-absorption is dom¬ 
inated by the strong transition involving the GT_ res¬ 
onance. With increasing A„, the cross section becomes 
insensitive to the energy dependence of the GT distribu¬ 
tion and depends only on the total GT_ strength. As a 
result we observe excellent agreement between the results 
of both approaches. 

Comparing the cross sections calculated with Pauli 
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blocking for the outgoing electron (see the lower panel 
of Fig. 9), we note that in [7], only allowed GT_ tran¬ 
sitions are taken into account when calculating neutrino 
absorption cross section. As is shown in Fig. 8, if the 
Pauli blocking is taken into account, forbidden down¬ 
ward transition become important with increasing tem¬ 
perature and electron chemical potential. It is apparent 
that such transitions along with the above discussed rea¬ 
sons lead to larger values for the TQRPA cross sections 
as compared to the shell-model results. 


IV. CONCLUSION 

In this work, we have studied thermal effects on the 
(anti)neutrino absorption for hot nuclei in supernova en¬ 
vironments. For this purpose, we have employed the 
proton-neutron QRPA extended to finite temperatures 
within the TFD formalism. As an example, cross sections 
were calculated for 56 Fe and 82 Ge in the temperature 
range from T = 0 to 1.72 MeV by taking into account the 
relevant charge-exchange transitions J 71 ' = , l =t , 2 =l= , 

and 3^. 

A detailed analysis of thermal effects was performed 
for allowed GT transitions which dominate the cross sec¬ 
tions for E u < 30 MeV neutrinos. Since the TQRPA 
does not support the Brink hypothesis, new peaks ap¬ 
pear in the GTip strength function at finite temperature 
due to transitions from the excited states. Moreover, 
thermal effects shift the GT resonance centroids to lower 
energies and this effect appears more strongly for the 
GT + strength in 82 Ge. The downward transitions from 
nuclear excited states were included in our calculations 
through detailed balance. The validity of detailed bal¬ 
ance for charge-exchange transitions is a consequence of 
the grand canonical treatment of hot nuclei. 

We have found that thermal effects on the GT strength 
enhance the absorption cross sections for low-energy 
(anti)neutrinos by several orders of magnitude. This en¬ 
hancement is mainly due to increasing contributions of 
downward transitions from excited states. However, in 
the supernova environment the electron chemical poten¬ 
tial increases more rapidly than temperature. As a re¬ 
sult, if the electron blocking in the final state is taken 
into account, the neutrino cross sections are drastically 
reduced. 

Although our calculations reveal the same thermal ef¬ 
fects as the shell-model calculation, the calculated low- 
energy cross sections for 56 Fe exceed the shell-model val¬ 
ues by two to three orders of magnitude. One of the possi¬ 
ble reason for this discrepancy is that the TQRPA under¬ 
estimates multinucleon correlations which are responsible 
for the GT strength fragmentation. On the other hand, 
the inclusion of back-resonances built on highly excited 
daughter states into shell-model calculations may also 
improve the agreement between the TQRPA and shell- 
model results. 

Since the TQRPA is not restricted by iron-group nu¬ 


clei, it has some advantages over shell-model calculations. 
To enhance its reliability and predictive power several 
improvements could be made. First of all, to account 
for multiconfigurational effects, the coupling of thermal 
charge-exchange phonons with more complex (e.g., two- 
phonon) configurations should be included into the ap¬ 
proach. At zero temperature this problem was consid¬ 
ered within the QPM [47, 48] by exploiting a separable 
form of the residual interaction. The other direction of 
the improvement is to combine the TQRPA method with 
self-consistent QRPA calculations based on either the rel¬ 
ativistic or Skyrrne nuclear energy density functionals. 
Recently, such calculations were performed at zero tem¬ 
perature [39, 41]. With a separable approximation for 
the Skyrrne interaction [49] it will be possible calculate 
the phonon coupling at T / 0 within a self-consistent 
theory. This is planned for the future. 
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Appendix 

Here we show that the detailed balance condition in 
the form (12) can be derived in a model independent 
way. When considering a grand canonical ensemble of 
hot nuclei, the probability to find the ith excited state of 
a nucleus with Z protons and N neutrons is given by 

P{£i,A%) = (2 Ji + 1) expj ——-—— }-Z _1 

(AT) 

where Z is the partition function of the grand canon¬ 
ical ensemble and J, is the angular momentum. No¬ 
tice that the excitation energies £-i are counted from the 
energy of noninteracting nucleons, i.e., Eq is a ground- 
state binding energy and the chemical potentials A„ jP 
do not include nucleon rest mass. Let us now introduce 
the temperature-dependent strength function for charge- 
exchange transitions as a thermal average of all transition 
strengths (probabilities) from states in the parent nucleus 
to states in the daughter nucleus: 

<F (T) ( E ) = 'Z'E P ^ A N)Slf ) S{E-Q\f). (A.2) 

Z,N ij 

Here Q^ = £f - £ t =F A M np and 

Slf = \(f,A z N % 1 1 \\T^\\i,A z N )\ 2 (A.3) 

is the reduced transition strength between states i and 
f in the parent and daughter nuclei, respectively. In the 
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above equations the upper sign corresponds to n —> p 
transitions, while the lower sign refers to p —> n transi¬ 
tions. The transition energy E can be both positive and 
negative. 

For the transition operators T^ and T^ + \ which dif¬ 
fer only by the isospin operator, the respective transition 
strengths, s\j and S^\ are connected by detailed bal¬ 
ance through 

(2J i + l)Slj ) = (2J f + l)S { + ) . (A.4) 

Combining this result with Eq. (A.2), we get the following 


relationship between the strength functions for n —> p 
and p —> n transitions in the thermal grand canonical 
ensemble 


$ (± )(-£) = exp{- E T (AA "P + AMnp) 

(A-5) 

This relation is exactly the same as derived within the 
TQRPA approach for charge-changing transitions in hot 
nuclei. 
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